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We present some exact results on bond percolation. We derive a relation that 
specifies the consequences for bond percolation quantities of replacing each bond 
of a lattice A by i bonds connecting the same adjacent vertices, thereby yielding 
the lattice A^. This relation is used to calculate the bond percolation threshold on 
A^. We show that this bond inflation leaves the universality class of the percolation 
transition invariant on a lattice of dimensionality d > 2 but changes it on a one- 
dimensional lattice and quasi-one-dimensional infinite-length strips. We also present 



analytic expressions for the average cluster number per vertex and correlation length 



1 for the bond percolation problem on the N — > oo limits of several families of N- vertex 

graphs. Finally, we explore the effect of bond vacancies on families of graphs with 
the property of bounded diameter as N — > oo. 

00 

o 

(N 



X 



I. INTRODUCTION 



Let G = G(V, E) be a connected graph defined by a set V of vertices ( = sites) and a 
set E of edges (= bonds) connecting pairs of vertices. Bond percolation is an interesting 
geometrical phenomenon in which one considers G to be modified in such a manner that 
each bond is independently randomly present with probability p. We denote the number 
of vertices and bonds as N = N(G) = \V\ and e(G) = \E\, and the formal limit of G as 
iV — y oo as {G} In the context of statistical mechanics, one is often interested in the 
iV — > oo limit of a regular lattice graph G, denoted A. As the bond occupation probability 
p decreases from 1, there are more and more absent bonds on A, and the probability P(p)a 
that a given vertex belongs to an infinite cluster decreases monotonically until, at a critical 
value, p c> Xi it vanishes and remains identically zero for < p < p C) &. This value, p Cj \, is 
the critical bond occupation probability, also called the critical bond percolation threshold. 
Other quantities also behave nonanalytically at p c ,h- For example, the average cluster size, 
S(p)a, increases monotonically as p increases from and diverges as p approaches p Cj A from 
below. Thus, the percolation transition is a geometrical transition from a region < p < p C) \ 
in which only finite connected clusters exist, to a region p c ^ < p < 1 in which there is a 
percolating cluster containing an infinite number of vertices and bonds. Another interesting 
quantity is the average number of (connected) clusters, including single vertices, divided 
by the total number of lattice vertices, in the limit N — > oo, denoted as (ti)a- Analogous 
statements hold for site percolation, in which each vertex of A is independently randomly 
present with probability p. Methods for studing percolation have included exact mappings, 
series expansions, Monte Carlo simulations, and the renormalization group. In addition to 
its intrinsic interest, percolation gives insight into a number of important phenomena such as 
the passage of fluids through porous media, electrical currents through composite materials 
consisting of conducting and insulating components, and the effect of lattice defects and 
disorder on thermal critical phenomena. Some reviews include [2]- (4] . 

A basic aspect of the percolation transition on a lattice A or, more generally, a limit 
{G} of a family of graphs, concerns the dependence of the critical behavior and the value of 
Pc{G} on properties of {G}. It is known that the percolation transition is in the universality 
class of the q = 1 ferromagnetic Potts model, and that its upper critical dimensionality 
is d u = 6. For a given lattice dimensionality d above the lower critical dimensionality, 
di — 1, the universality class of the percolation transition is independent of details of the 
lattice structure, such as the coordination number. In contrast, other quantities, such as the 
critical bond percolation value p c ,A, are non- universal and do depend on lattice properties like 
coordination number The percolation threshold p C: \ is usually a decreasing function 
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of the lattice coordination number, although certain exceptions have been found 10l. 111]. 

In this paper we report some exact analytic results concerning the factors that deter- 
mine the critical behavior of the percolation transition and the dependence of percolation 
quantities on the properties of {G}. Specifically, we derive a relation that specifies the con- 
sequences for bond percolation quantities of replacing each bond of a lattice A by £ bonds 
connecting the same adjacent vertices. This result elucidates the effect of an arbitrarily 
great increase in the degrees of the vertices of G and its N — > oo limit, {G}. (Here, the 
degree of a vertex is defined as the number of bonds connecting to this vertex.) We show 
that this £-fold bond inflation leaves the universality class of the percolation transition in- 
variant on a lattice of dimensionality d > 2 but changes it on a one-dimensional lattice and 
on quasi-one-dimensional infinite-length, finite-width lattice strips. This is demonstrated by 
changes in the critical exponents governing the divergences in the correlation length and 
average cluster size as p /* 1. We also present analytic expressions for the average cluster 
number per vertex, (n) , for some families of graphs containing repeated complete graphs K r 



151 ] as subgraphs. As part of this, we discuss the analytic properties of (n) in the complex 



p plane and use these to determine the radius of convergence of the relevant small-p series 
expansion of (n) . Finally, we explore the effects of bond vacancies on families of graphs with 
the property of bounded diameter as N — y oo. The present work extends our previous 
study in [17 1 . 



II. GENERAL BACKGROUND AND CALCULATIONAL METHODS 

In this section we discuss the calculational methods that we employ. We make use of the 
fact that several quantities in bond percolation can be obtained from the partition function 
for the g-state Potts model 0, 18-20 1 in the limit q — >■ 1. We review this connection next. 



In thermal equilibrium at temperature T, the general Potts model partition function in an 
external magnetic field H is given by Z = ^\ , with the Hamiltonian 

(ij) I 

where i, j, I label vertices of G, G\ are classical spin variables on these vertices taking values 
in the set I q = {1, ...,(?}, = (/c^T) -1 , and (ij) denote pairs of adjacent vertices. 

The zero-field Potts model Hamiltonian % and partition function Z are invariant under 
the global transformation in which <jj — > gcXi V i G V, with g G S q , where S q is the symmetric 
(= permutation) group on q objects. Because of this invariance, one can, without loss of 
generality, choose the value picked out by the external magnetic field H to be = 1, as we 
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have done in ( 12.11) . It will be convenient to introduce the notation 

K = (3 J , h = (3H, a = e K , v = a-l, w = e h . (2.2) 

Given a graph G = (V,E), a spanning subgraph G' C G is defined as the subgraph 
containing the same set of vertices V and a subset of the bonds of G; G' = (V, E') with 
E 1 C i£. We denote the number of connected components of G and G' as n(G) and n(G'), 
respectively. The property that G is connected is the statement that n(G) = 1. The link 
between percolation and the Potts model stems from the property that the Potts model 
partition function can be expressed as a purely graph-theoretic sum over contributions from 



spanning subgraphs G 1 . For H = 0, this expression is [18] 

Z(G,q,v)= 5> £ ( G V (G ' } . (2.3) 

G'CG 

As is evident from (12. 3p . Z(G,q,v) is a polynomial in q and v. This expression also allows 
one to generalize q from positive integers to real numbers. 

As an example of the calculation of percolation quantities from (12. 3p . we consider the 
average number of clusters per vertex. On a graph G, this is given by 



<n> 



Y. G >v e{G '\i-v) 



e(G)-e(G>) 



(l/N)j: G ,n(G')[p/(l-p)r^ 

£ G <b/(i-p)] e(G,) ' l ' J 

This follows because each G' contains n(G') connected components, and appears in the 
numerator of the expression in the first line with weight given by p e ^ G '\l — p) e ( G )~ e ( G ") ) since 
the probability that all of the bonds in G' are present is p e ( G "> and the probability that all of 
the other e(G) — e(G') bonds in G are absent is (1 -p)< G )~< G '). This sum in the numerator 
over the set of spanning subgraphs G' is normalized by the indicated denominator and by 
an overall factor of 1/7V to obtain the average number of connected components (clusters) 
per vertex. 

On a finite graph G one defines the (reduced) free energy per vertex of the Potts model 

as 

f(G,q,v)=\n[Z(G,q,v)^ N } (2.5) 

and, in the limit iV — > oo, 

f({G},q,v)= lim f(G,q,v) , (2.6) 

jV->oo 
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If one sets 



v = v p = - , i.e., p = —— = 1 - e , (2.7) 

1 — p v + 1 

differentiates f(G, q, v p ) with respect to q, and then sets q — 1, one obtains (n) n , as given in 

eq. (El, i.e., 

df(G,q,v p ) 



(n) 



N 



In particular, in the limit iV — >■ oo, 



(n> 



df({G},g,v p 



<?=i 



<9g 



(2.8) 



(2.9) 



The formula (I2.9p relates a geometric property of (the N — > oo limit of) a bond-diluted 
graph with the derivative of the reduced free energy of the zero-field Potts model, evaluated 
at a certain temperature, as q — > 1, on a graph with no bond dilution. The mapping 
(12. 7p . in conjunction with eq. (I2.9p . formally associates the interval < p < 1 with the 
interval < v < oo, which is the physical range of values of the temperature variable v for 
the ferromagnetic Potts model, with p — > + corresponding to temperature T — > oo, i.e., 
v + , and p — > 1 to T — > 0, i.e, t> — > oo. 

For the case of nonzero field, denote the connected subgraphs of a spanning subgraph 
G' as i = l,..,n(G"). Then one can obtain a generalized expression for the partition 
function as a sum over contributions from spanning subgraphs, as 19] 



n(G') 



Z(G,q,v,w)= II^- 1 ^^) 

G'CG i=l 



(2.10) 



Some general properties of Z(G, q, v, w) were derived and exact results for families of graphs 



2l|-[25|. For the zero-field special case, this reduces to Z(G, q,v, 1) = 



were given m 
Z(G,q,v). 

Extending the zero-field definition, we define the dimensionless reduced free energy of the 
Potts model in an external field as 



f({G},q,v,w) 



lim —]n[Z(G,q,v,w)] 



(2-11) 



For a given {G}, the quantities P(p) and S(p) can then be determined in terms of deriva- 
tives of this reduced free energy: 



P{p) 



1 + 



d 2 f 
dhdq 



5=1, h=0+ 



l + W 



d 2 f 

dwdq 



q=l, w=l~' 
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(2.12) 



and 



S(p) 



d 3 f 



(dh) 2 dq 



=1, h=0+ 



= w JL\ w JL(?lX\ ( 2 13) 

dw L dw \dqJl q=i, w=i+ 

where h = + and w = 1 + mean lini/j^ + an d lim w _ + x+- 

We recall the forms of the singularities in some quantities at the percolation transition. 
On a given lattice A, as the bond occupation probability increases from to 1, (n)\ decreases 
monotonically from 1 to and has a singularity at the critical threshold value p = p c ^ of 
the form 

(n)A,sing. ~ \p -Pc,\\ 2 ~ a as p->j) c ,A . (2.14) 
As p decreases from 1 to p c ,A, the average probability that a vertex is connected to the infinite 
cluster, P(p)a, decreases to zero like 

P{p) A oc (p - p cA f as p \ p cA ■ (2.15) 

As p increases from toward p c ,A, the average size of a percolation cluster, S(p)\, diverges 
like l/(p c ,A — p) r ' ■ As p decreases toward p Cj A from above, S(p)a, defined as the average 
size of finite clusters, diverges like l/(p — Pca) 1 ■ Scaling and renormalization group methods 
yield 7 = 7', so that 

S(p)a oc as p -> p C: A ■ (2.16) 

\p-Pc,Ar 

Similarly, one can define a correlation length £ (p) characterizing the size of clusters, and this 
diverges at the percolation transition like 

£(p) « 1 — vz as P->Pc,a ■ (2.17) 

In the Potts model, the divergence in the correlation length corresponds to the equality of 
the leading and subleading eigenvalues of the transfer matrix, 

* = Vn 7x V • ( 2 ' 18 ) 

HlyA max / ^submax ) 

Setting q — 1 and v — v p in Eq. (12. 7B . one thus determines the divergence in £(p) for the 
bond percolation problem as p — > p Ct A from the divergence in £ for the ferromagnetic Potts 
model. 
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The relation of percolation to the q = 1 Potts ferromagnet implies, in particular, that 
the percolation transition on two-dimensional lattices is in the same universality class as the 
two-dimensional q = 1 Potts ferromagnet, with the exactly known exponents a = —2/3, 
(3 = 5/36, 7 = 43/18, v = 4/3, etc. [20|. In the context of conformal algebra, these critical 
exponents are associated with a rational conformal field theory with central charge c = 



(e.g., |26|). The percolation transition on a (^-dimensional lattice with d ^ 2 is in a different 
universality class. In particular, as noted before, the upper critical dimensionality for the 
percolation transition is d u = 6, corresponding to a Ginzburg-Landau function with highest 
power 3 and critical exponents a = —1, (3 = 1, 7 = 1, and v = 1/2. This relation is 
embodied in the exact solution on the Cayley tree (reviewed, e.g., in |2J). 

We denote the critical point in K of the g-state Potts ferromagnet on the lattice A as K Ctq ^, 
and the corresponding value of v as f c , ? ,A- We also denote -^ Ci i,A = hnig-^i K c , q ,A- The critical 
percolation occupation probability, i.e., the percolation threshold, p Ci A, is determined in terms 
of the critical point of the corresponding q — > 1 Potts model as p c .A = 1 — e~ Kc ' 1 - A . Exact 
values of p c A for bond percolation on the square (sq), triangular (tri), and honeycomb (he) 
lattices were obtained in Q| (reviews include [^-{ij]). For later reference, these well-known 
critical bond percolation thresholds are p CtSg = 1/2, p Cjtr i = 2sin(7r/18) = 0.347296.., and 
Pc,hc — 1 —Pc,tri — 1 — 2sin(7r/18) = 0.652704.. Much work has been done on determinations 
of p c ,A values for various lattices 0]- 14]. For example, for the simple cubic (sc) lattice, 
Pc,sc = 0.2488126(5) [7j], where the number in parentheses is the uncertainty in the last digit. 
After reviewing this background, we now proceed to our new results. 



III. BOND PERCOLATION QUANTITIES ON BOND-INFLATED LATTICES 

A. Basic Result 

Consider an arbitrary graph G = (V, E). We define the graph Ge as the graph obtained 
by replacing each bond of G by i bonds connecting the same vertices. Similarly, in the 
iV — > 00 limit, we define {G}i in the same manner. A K-regular graph G is defined as a 
graph all of whose vertices have the same degree (coordination number), k. To cover the 
case of graphs that are not K-regular, it will also be useful to define an effective vertex degree 
in the — > 00 limit, namely 

«e//({G}) = hm , (3.1) 

For example, this can be defined on duals of Archimedean lattices j^]. Clearly, for a k- 
regular graph G, n e ff = k. For a k- regular graph G, the Mold bond inflation increases the 
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vertex degree to I k, and similarly, for graphs that have vertices of different degrees, the 
degree of each vertex is increased by the factor i. These cases are subsumed in the N — > oo 
limit, as 

Keff({G})t = tKeff{{G}) . (3.2) 

In particular, if {G} is a regular lattice A, we denote the lattice Ai as the result of replacing 
each bond on A by t bonds connecting the same adjacent vertices. 

In this section we derive exact relations that express quantities characterizing bond per- 
colation on {G}g in terms of the corresponding quantities on the original {G} with a trans- 
formed value of the bond occupation probability. The starting point of our derivation is the 
observation that the effect of £-fold bond inflation is embodied in the relation 

Z(Gi, q, v) = Z(G, q, ve) , (3.3) 

where 

v t = (v + 1) £ - 1 . (3.4) 



Some other implications of this are discussed in [28| . Because the interaction of the external 
magnetic field with the spins crj in (12. ip is unaffected by the bond inflation, the generalization 
of Eq. (13. 3p to the case H ^ is immediate: 

Z(G e , q, v, w) = Z(G, q, v £ , w) . (3.5) 

For our particular application, we set v = v p in Eq. (12 .7p . which yields the relation 



p, = l-(l_p)<=p W'V_p)i-i 

7=1 



(3.6) 



where (.) = £\/[j\(£ — is the binomial coefficient. The first few explicit cases, aside from 
Pi = p, are 

p 2 = p(2-p), (3.7) 

p 3 =p(3-3p + p 2 ) , (3.8) 

Pi = p(4 — 6p + Ap 2 — p 3 ) 

= p(2-p)(2-2p + p 2 ) , (3.9) 

and so forth for higher I. 
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From Eq. ( 13. 6j) it follows that the r'th derivative of pt with respect to p is 

-iy-H (r) (i - P y- r (3.10) 



for 1 < r < £ and zero for r > £ + 1, where £m is the falling factorial, 

j (r) = J](J- a ) . (3.11) 

We will need to invert Eq. (13.61) and solve for p in terms of pe. For this purpose, we note 
that Eq. (13. 6p is an tth degree algebraic equation for p, but the relevant one among the £ 
roots is determined by the requirement that v — if and only if ve = 0, so p = if and only 
if pi — 0. This root is given by 

p = l-(l-p i ) 1 ". (3.12) 

Note that this relation has the same form as Eq. (13.61) with the replacements p -H- pe and 
£ ^ l/£. 

From Eqs. (I3.6P and (13. 12[) we derive the following properties. First, the transforma- 
tion (13.61) maps the interval p G [0, 1] to the interval pe G [0, 1], and similarly the inverse 
transformation (13.121) maps the interval pi G [0, 1] to p G [0, 1]. Second, 

p = <{=► p e = (3.13) 

and 

p = \ ^ Vl = \ . (3.14) 

Third, for £ > 2, 

^ - p = (1 -p)pe~i (3.15) 
Fourth, as is evident from Eq. (I3.15p . 

Pt>p for p G [0, 1] and £ > 2 , (3.16) 

with equality only at p = pi = and p = pe = 1. (For £ = 1, Pl — p, so this inequality 
is realized as an equality for all p.) For £ > 2, the difference pe — p has a maximum in the 
interval < p < 1 which starts at p = 1/2 for £ = 2 and moves to the left as £ increases. 
Fifth, for fixed p6 (0,1), 

Pi is a monotonically increasing function of £ for p G (0, 1) (3-17) 

and for fixed P £ G (0, 1), 

p is a monotonically decreasing function of £ for pi G (0, 1) . (3.18) 
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Sixth, from the r = 1 special case of (13.10p . dpi/dp = £{l — p) 1 " 1 , it follows that for fixed £, 
Pi is a monotonically increasing function of p for p e (0, 1) . (3.19) 

Similarly, for fixed £, 

p is a monotonically increasing function of pe for p e (0, 1) . (3.20) 

For £ > 3, the curve of pi as a function of p is quite flat near p = 1 because, as is evident 
from Eq. (13.101) . the second derivative d 2 pe/dp 2 vanishes at p = 1. More generally, the r'th 
derivative, d r pi/dp r , vanishes at p = 1 for £ > r + 1. 

From Eqs. (13. 5ft and (13. 6ft . it follows that an arbitrary percolation quantity Q{j>){G} e on 
{G}i, such as {n){G} e , Pij>){G}^ S(p){c} e , etc. satisfies the relation 

Q(p){G h = Q(Pi){G} ■ (3.21) 

In particular, this relation holds for iV — > oo limits of lattice graphs {G} = A. Thus, for 
example, 

(n)(p)A e = (n)( Pe ) A , (3.22) 

P(p) Ae = P( Pe ) A , (3.23) 

S(ph e = S(p e ) A , (3.24) 

C(p)a, = £(P/)a , (3.25) 

and so forth for other bond percolation quantities on the bond-inflated lattice A^, where pi 
is given in terms of p by Eq. (13.61) . From either (I3.23P or (I3.24p in conjunction with (13. 6p . it 
follows that 

p c ,A e = I - (I - Pc,A) 1/e (3.26) 

and hence, in particular, 

Pc,A e < Pc,A , (3.27) 

with equality only for the case where p CtA = p c ,A e = 1- These relations (13.61) . (I3.17p . (I3.2ip . 
and (13.261) are important exact results, since they describe the effect of the Mold bond 
inflation on percolation quantities. This £-fold bond inflation increases the vertex degree by 
the factor £ and has the consequence that percolation quantities on {G}i are equal to the 
corresponding quantities on {G} with the bond occupation probability p replaced by the 
larger probabiity pi. This means that as p increases, the infinite percolation cluster appears 
at a smaller value of p on {G}e than on {G}, as given by Eq. (13.261) . That is, the Mold bond 
inflation enhances the formation of an infinite percolating cluster on the resultant graph 
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{G}e. The relation (13.261) is specific to this process of £-fo\d bond inflation, while the general 
inequality (I3.27P follows from the fact that if a lattice A' is obtained from a lattice A by the 
addition of (an arbitrary) set of bonds, then p c ^> < p c ,A- Recall that this fact is clear, since 
if p is large enough for an infinite percolation cluster to exist on A, then this cluster certainly 
also exists on A' with its additional bonds. 

The critical percolation threshold depends not only on the lattice dimension d and vertex 
degree k for a regular K-regular lattice A or effective vertex degree n e ff for a lattice with 
vertices of different degrees, but also on other related properties of the lattice. Several studies 
of lattice properties that affect p c A have been carried out for regular lattices, including 
Archimedean lattices and their planar duals Q-jl^- Clearly, another graphical property 
that is relevant is the edge-connectivity (= bond-connectivity) A(G), defined as the minimum 
number of bonds that must be removed to increase the number of components by one. In 
particular, given that G is connected, i.e., n(G) = 1, the edge connectivity, A(G) is the 
minimum number of bonds that must be removed to separate the graph into two disconnected 
components. This is related to, but different from, vertex degree measures such as k for a 
K-regular graph or K e ff for the limit N — > oo of a graph with vertices of several different 
degrees. For example, a tree graph G t (defined as a connected graph with no circuits) may 
have vertices with arbitrary degrees, but has X(G t ) = 1. For a K-regular lattice with periodic 
boundary conditions, A(A) = k. For our present discussion, we note that 

\(G £ )=e\{G). (3.28) 

Thus, for an arbitrary graph G, the £-fold bond inflation increases both the vertex degree 
and the edge-connectivity by the factor of I. 

Although only the behavior of p and pi in the interval pG [0, 1] and thus pi G [0, 1] is of 
direct interest for percolation, the behavior for real p and pe outside this interval is also of 
interest in a broader context. We note that for p < 0, pe is also negative and for p » 1, pe 
is positive if I is odd and negative if I is even. Furthermore, for I > 2, the point p = 1 is a 
(global) maximum for pi if I is even and an inflection point if i is odd. Moreover, in addition 
to its zero at p = 0, pe vanishes at p = 2 if and only if £ is even. Even more generally, as will 
be discussed below, the analytic behavior of various percolation quantities in the complex p 
plane is of interest. Indeed, there are cases where the radius of convergence of a Taylor series 
expansion for (n) about the point p = is not set by p c , but instead by complex singularities 



in the p plane. Explicit examples of this were given in Ref. 17] . 
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TABLE I: Values of the critical bond percolation threshold probability p c .A e on the lattice Ki obtained 
from the lattice A by replacing each bond by I bonds (so that Ai = A). We list results for the square (sq), 
triangular (tri), honeycomb (he), and simple cubic (sc) lattices. See text for further details. 



ft 


Pc,(sc) e 


Pc,(tri) e 


Pc,(sq) e 


Pc,(hc) e 


1 
1 




n "3/1 7 


u.ouu 


U.OOo 


2 


0.133 


0.192 


0.293 


0.411 


3 


0.0910 


0.133 


0.206 


0.297 


4 


0.0690 


0.101 


0.159 


0.232 


5 


0.0556 


0.0818 


0.129 


0.191 


6 


0.0466 


0.0686 


0.109 


0.162 


7 


0.0400 


0.0591 


0.0943 


0.140 


8 


0.0351 


0.0519 


0.0830 


0.124 



B. Percolation Threshold on Some Specific Bond-Inflated Lattices 



We now focus on regular lattice graphs {G} = A and {Ge} = Ae. It is of interest to apply 
our general result f!3.26j) to obtain some illustrative numerical values of p c ,K t on various 
lattices. First, on a one-dimensional or quasi-one-dimensional (Q1D) infinite-length, finite- 
width strip Aqid, with j> c ,qid = 1, the Mold bond inflation leaves this property unchanged, 
i.e., Pc,A Q1D = 1, as is clear from ( 13. 6 p and ( I3.14p . On higher- dimensional lattices, as special 
cases of our general relation (13. 26 p . we display the following illustrative results: 



Pc,(sq)i 



(1 - Pc,sq) 



In 1 
2, 



Pc,(tri) e 



Pc,(hc) e 



(1 - Petri) 



1 - 2 sin 



i/t 



1 - (1 - p cM ) 
1-fl 



2 sin 



7T 



l/t 



(3.29) 

(3.30) 
(3.31) 



p cAsc)l = l-(l- p CjSC p~ = l - (0.751187) 1/£ . (3.32) 

We list numerical values of these threshold percolation probabilities in Table [H Note that 
n((hc)e) = 3£, K((sq)e) = 4£, and K((tri)e) = k((sc)^) = 6£. 
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C. Effect of Bond Inflation on Universality Class of Percolation Transition 



Consider a regular d- dimensional lattice A with d > 2, so that p Cj A € (0,1). As noted 
above, the percolation transition is in the same universality class as the phase transition 
in the ferromagnetic Potts model in the limit q — > 1. The universality class of a finite- 
temperature ferromagnetic phase transition depends only on the symmetry group of the 
(zero-field) Hamiltonian and the dimensionality of the lattice. In particular, this univer- 
sality class is independent of the coordination number of the lattice. From these facts it 
follows that the process of bond inflation does not change the universality class of the per- 
colation transition on such a lattice. In contrast, special properties apply in the case of a 
one- dimensional lattice and infinite-length, finite-width lattice strips, which are quasi-one- 
dimensional. For these there is no finite-temperature phase transition in a spin model (with 
short-range spin-spin interactions), and, correspondingly, the critical percolation threshold 
is p c — 1. We will show below that for one-dimensional and quasi-one- dimensional lattices, 
bond inflation can change a critical exponent characterizing percolation. 



D. Effect of Bond Inflation on Percolation on a ID Lattice 

We denote the ID lattices of length A" vertices and free and cyclic boundary conditions 
as L n and C n , respectively and the corresponding Mold bond-inflated lattices as (L n )g and 
(C n )g. The A^ — > oo limits of these lattices are denoted {L}, {C}, {L}e, and {C}i. Since 
percolation quantities are independent of the boundary conditions we will denote both of 
these limits simply as ID, and the corresponding £-fold bond-inflated limit as (\D)g. It 
will be convenient to recall how the well-known results for percolation on a one-dimensional 
lattice follow from the solution to the Potts model in the q — > 1 limit. Since this model is 
only critical at T = 0, i.e., v = oo, it follows via Eq. ( 12. 7p that 

Pc,id = 1 , (3.33) 

Evaluating Eq. ( 12. 9 j) with the reduced free energy f 1D = ln(g + v), one has 

(n) 1D = l-p, (3.34) 

For any p G [0, 1), in the limit A" — > oo, there is zero probability that a vertex is in an infinite 
cluster because there is no infinite cluster. Such a cluster only exists for p = p c> \o = 1- Hence, 

ifpe[0,l) 

P{v)id={ , (3.35) 

1 if p = 1 
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This is analogous to the singularity in the magnetization for a ID spin model, which is 
identically zero for any finite temperature and jumps to 1 at the critical temperature, T = 0. 
The ID Potts model correlation function has the form 

G{r) oc p r , (3.36) 

where p is the ratio of the next-to-maximal eigenvalue of the transfer matrix to the maximal 
eigenvalue, 

^■submax /o nrr\ 

p = — . (3.37) 

Here, pm = v/(q + v). Setting q = I and v = v p in Eq. (12. 7p gives p = p, so 

G(r) 1D =p r , (3.38) 
With G(r) ~ e~ r ^ for r — > oo and p ^ p c , one has 

$id = -7— (3.39) 
mp 

and hence 

Zid ~ 7^ as p / 1 . (3.40) 
1 — p 

Therefore, the corresponding critical exponent is 

^id = 1 ■ (3.41) 

Similarly, 

S(p)iD = ^ , (3.42) 
1 — p 

so that S(p)iD diverges as p /* 1 with the critical exponent 

7id = 1 ■ (3.43) 

Having reviewed these well-known results, we now analyze the effects of £-fold bond 
inflation. From Eqs. ( 13. 6 j) or ( I3.26p . it follows that 

Pc(iD), = 1 • (3.44) 

As special cases of our general result (|3.21|) . we have 

(n) (1D)e = (l- P y , (3.45) 
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ifpe [0,1) 

P{p\ 1D)t = P(p) {m = t , (3.46) 

1 it p — 1 

and 

As p /* p Cj i_D — 1) this diverges like S(p) ~ 2/(1 — p) . From this we find that the critical 
exponent for bond percolation on the Mold bond-inflated lattice (lD)e is 



Furthermore, 



so 



Therefore, 



7 (1D) , = I • (3.48) 
G(r) (1D)e = (p e y , (3.49) 
«"* - -5^ (3-50) 

^~i^ = (T3^' asp/1 (3 - 51) 

and hence the correlation length diverges as p / 1 on the Mold bond-inflated lattice (1-D)^ 
with critical exponent 

V{iD) t = I ■ (3-52) 

These are important results, because they show that, in contrast to percolation on higher- 
dimensional lattices, where p c G (0, 1), here the bond-inflation changes the critical exponents 
and hence the universality class of the percolation transition. 



E. Effect of Bond Inflation on a Quasi-One-Dimensional Lattice Strips 

Our result that Mold bond inflation changes the universality class of the percolation 
transition on a one-dimensional lattice generalizes to apply also to quasi-one- dimensional, 
infinite-length, finite-width strips {G s }. This is a consequence of the fact that, because the 
Potts ferromagnet is only critical at T = (i.e., v = oo) on such strips, corresponding to 
Pc,{G a } = 1) divergences of the form 

and 

«")«« <* (1 -,)-.-. < X54 > 
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as p 1 change to 

S(p){G ^ K (i-pV {G -> = (i- P W> (3 - 55) 

and 

e(p){G ->< K (i- P ]ync s} = ( i-p)^ G .> • (3 - 56) 

Hence, the corresponding critical exponents on the £-fold bond-inflated infinite-length, finite- 
width strip graph {G s }e are 

l{G a } e = h{G a } (3.57) 

and 

V{G s } e = Zv{G s } • (3.58) 

Thus, again the universality class is changed by the £-fold bond inflation on these infinite- 
length, finite-width strip graphs. 



IV. PERCOLATION ON SPECIFIC INFINITE-LENGTH LATTICE STRIP 

GRAPHS 

In this section we give some results for the infinite-length limits of lattice strip graphs. 
We begin with the L x — > oo limit of the square-lattice strip graph with width L y = 2 and 
free (F) transverse boundary conditions, i.e., the ladder graph. The longitudinal boundary 
conditions are free or periodic. We denote this limit as sq,2p, where the subscript F refers 
to the transverse boundary conditions. We have previously calculated the average per-site 



cluster number, which is [17] 



;i -p) 2 (2 + p-2p 2 ) 
2(1 -p 2 + p 3 ) 



\ n )sq,2 F = —i Z3 , ^ • (4- 1 ) 



A notable property of this exact result is that it has a pole singularity at a negative real 
value of p, namely p = —0.7549 (as well as two complex values of p) and, as is evident, 
the singularity at this negative real value is closer to the origin than the critical percolation 
threshold value, p c — 1. Hence, the radius of convergence of a small-p Taylor series expansion 
of this quantity is set not by p c , but by this unphysical singularity. This phenomenon of an 
unphysical singularity being closer to the origin than the physical singularity is also true of 
our result in Eq. (I5.14j) below and of many infinite-length, finite-width strips analyzed in 



13]. 
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Here we go on to calculate the divergences in £(p) and S(p) for this strip as p — y p c = 1. 
Using the solution in Ref. 29] for the general Potts model partition function on a strip of 
this type of arbitrary length, we calculate 

^sq,2p ,submax t A n\ 

Psq,2 F = — T , 

^sq,2p ,max 

where 



A 



sq,2p,max 



v 3 +4v 2 +3qv+q 2 + 



v +Av°-2qv' i -2q 2 v 3 +l2v' i +l6qv d +l3q z v z +6q- i v+q' 1 



and 



A 



sq,2p,submax 



q + v + 4) + 



v 4 + 4w 3 + 12f 2 - 2gv 2 + Aqv + g 2 



lV2l 



In the limit K — >■ oo, i.e., a — )■ oo, this leads to a divergence in £ of the form 



£so,2 P ~ — as a — > oo 



1/2 

(4.3) 
(4.4) 

(4.5) 



where a = e K , as defined in (12. 2p . Setting g = 1 and v — v p as in Eq. (12 .7p . we thus obtain 

1 



e 



sq,2p 



[i-pY 



as p / 1 



(4.6) 



Hence, the corresponding correlation-length critical exponent for bond percolation on the 
infinite ladder graph is 



(4.7) 



By similar methods we derive the following results for some other infinite-length limits of 
finite- width lattice strips. As before, our results apply for either free or periodic longitudinal 
boundary conditions. First, we consider the strip of the square lattice with (transverse) 
width L y = 2, with periodic (P), rather than free, transverse boundary conditions. In effect, 
this doubles each transverse bond, leaving the longitudinal bonds unchanged. We denote 
the — y oo limit of this strip as sq,2 P . (n) sq ^ P was calculated in [17|. Using our results 



from 



30|, we calculate the divergence in the correlation length for the ferromagnetic Potts 



model as T — y on this strip to be 



£sg,2 P ~ — as ci y oo . 
' ' q 



(4- 
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Setting q = 1 and v = v p , we derive the corresponding results for the divergence in £ for 
bond percolation on this strip: 

Zaq,2p OC ^ j ^ 2 aS p / 1 , (4.9) 

which is the same as for the sq,2p strip. Hence, the correlation- length critical exponent for 
bond percolation on this strip is 

u sq , 2p = 2 . (4.10) 
Next, we consider the width L y = 2 strip of the triangular lattice with free transverse 



boundary conditions, denoted tri.2 F . We calculated (n) tr i,2 F hi [17]. Here, using our results 
in j3l|, we calculate the divergence in the correlation length for the Potts ferromagnet on 
this strip as T — > to be 

a 3 

im,2 F ~ 7T as a oo (4.11) 
Zq 

and hence for the bond percolation problem, 

Wf oc as p / 1 . (4.12) 

Consequently, 

^«,2 F = 3 . (4.13) 

These results again show how the critical behavior of percolation is sensitive to details 
of the lattice structure at the lower critical dimensionality. Note that Km = 2, K S q,2 F — 3, 
K S q,2 P = 4, and Ktri,2 F = 4 for these strips with periodic longitudinal boundary conditions. 
(More generally, for the corresponding strips with free longitudinal boundary conditions, 
these values apply for the respective K e /f.) Carrying out Mold bond inflation on these 
infinite-length strips, one again gets a change in the critical exponents and hence universality 
class for the percolation transition, as a special case of (I3.58p . 



V. BOND PERCOLATION ON G[K r ,jn] 
A. General 

In this section we present exact expressions for the average cluster number (n) for 
the infinite-length limits of a family of graphs with variable vertex degree, namely 
lirn m _j.oo G[(K r ) m , jn, BC], defined as the infinite-length limit of a line or ring of m sub- 
graphs K r connected in such a manner that all vertices of the tth K r are connected to 
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all vertices of the (i + l)'th K r . Here the notation BC refers to the longitudinal bound- 
ary conditions, which are free (FBC) for the line and periodic (PBC) for the ring. We 
also present results for the divergence in the correlation length. One of the reasons that 
this family of graphs is useful for the present study is that one can vary the vertex degree 
over a rather wide range by varying r. An illustrative example of a member of the family 
G[(K r ) m , jn, PBC] for the case r = 2 and m = 4 was given as Fig. 1(a) in Ref. 32). The 
cyclic strip G[(K r ) m , jn, PBC] is a ^-regular graph with uniform vertex degree 

k = 3r — 1 for G[(K r ) m ,jn,PBC\ , (5.1) 

and all of the vertices except the end vertices of the free strip G[(K r ) m , jn, FBC] also have 
this degree. Thus G[(K r ) m , jn, PBC] is a (3r — l)-connected graph. For both FBC and 
PBC, the graph G[(K r ) m , jn, BC] has iV = mr vertices so that the infinite-length limit can 
be written equivalently as m — > oo or N — > oo. Our result for (n) depends on r, but not on 
these boundary conditions. We denote the formal limit m — > oo of this family as 

G[K r ,jn] = lim G[(K r ) m ,jn] (5.2) 



where we suppress the BC in the notation, since (n) is independent of the boundary condi- 
tions. For FBC, as well as PBC, in the m — > oo limit, we have 

K e jf = 3r — 1 for G[K r ,jn] . (5.3) 



Our method to calculate (n) in [17J and here is to apply eq. (12. 9p in conjunction with 
exact results that we have computed for the free energy of the Potts model on infinite- 
length, finite-width strips of various lattices. We will give explicit expressions for the cases 
r = 1,2,3 and relevant properties of the result for general r. We will also display Taylor 
series expansions of the resultant (n) for p near and for p near 1, in the latter case, using the 
expansion variable s = 1 — p. Our current work is a continuation of our previous calculations 



of (n) for various families of graphs in Ref. [17|. Other studies of (n) for bond percolation 



include 



33]. 



B. Calculations for G[K r ,jn] 

For the calculation of (n), one needs the reduced free energy for the m — > oo limit of 
the strip graph G[(K r ) m , jn, BC]. This is simplest for the case of free boundary conditions. 
In 34j we determined the general structural form of the Potts model partition function 



Z(G[(K r ) m ,jn,BC],q,v) for this family. As we discussed in [17|, to calculate (n) on the 
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iV — > oo limit of a one-parameter family of recursive graphs, such as lattice strip graphs, using 
the free energy of the Potts model on these graphs, one needs the dominant term contributing 
to the free energy in the ferromagnetic region. For these classes of one-parameter graphs 
G m of length m subunits, the free energy has the form of a sum of m'th powers of certain 
functions, multiplied by certain coefficients of degree d, ranging from to a maximal degree d 
depending on the transverse structure of the strip but not its length. The dominant term in 
the ferromagnetic region arises from the degree d = sector, and is the same independent of 
the longitudinal boundary condition, in accord with the requirement that the thermodynamic 
behavior should be independent of the boundary conditions in the infinite-length limit. For 
;he present case, we proved that Z(G[(K r ) m , jn, PBC], q, v) has the general structural form 



34] 



nx(r,d) 

Z(G[(K r ) m ,jn},q,v) = J2^ d Yl ( X G[K r ,MAi) m , (5-4) 

d=0 i=l 

where 

r-d+l / 1 

r — 1 



n T (r,d)= Li) ( 5 - 5 ) 

and the coefficient fid is a polynomial in q of degree d given by 



fio = 1 (5.6) 

and 

where qu\ is the falling factorial defined in (13. lip . (The symbols nx{r,d) in Eqs. (15 .4p and 
( 15. 5 p above and nzh{.L y ,Gr,,d) in Eq. ( 18. ip follow the notation used in our earlier papers 
and should not be confused with the notation n(G') for cluster numbers.) 
With FBC, only the \G[K r ,jn],o,i contribute. For general r, we have 

f{G[K r ,jn],q,v) = - ki[\ G[Kr!jn]Amax \ , (5.8) 

where \G{K ri jn],o,max denotes the \G[K r jn],o,j of maximal magnitude for ferromagnetic v = v p , 
i.e., real positive v corresponding to p G [0, 1]. Thus, 

df(G[K r ,jn],q,v p ) 
(n) G [K r ,jn] = — . • (5-9) 



dq 



q=l 



If r = 1, the graph G[(Ki) m , jn, BC] is the path graph L m on m vertices for free boundary 
conditions, and the circuit graph C m for periodic BC, which has already been discussed above. 
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C. G[K 2 ,jn] 



Here we calculate the average cluster number, per vertex and the correlation length for the 
m — > oo limit of the family G[(K 2 ) m , jn, BC], denoted as G[K 2 ,jn]. A graph in this family- 
can also be regarded as a square-lattice ladder graph with next-nearest-neighbor bonds. 
From our calculation of the partition function for this graph in 35j, we have 



f{G[K 2 ,jn],q,v) = ^\n[X G[K2:jn]:+ ) 



(5.10) 



where 



with 



A 



G[K 2 ,jn},± 



T K2 ± ^/R 



T, 



K 2 



v 5 + 5w 4 + 10v 3 + 12v 2 + 5qv + q 2 



(5.11) 



(5.12) 



and 



R K2 = v w + I0v 9 + 45v 8 + 116f 7 + 196v 6 + 224t> 5 + lUv 4 - 6v 6 q - 10v 5 q + 40t> 4 g + 104v 3 q 

- 2v 5 q 2 - 10v 4 q 2 - Av 3 q 2 + Alv 2 q 2 + IQvq 3 + q 4 . (5.13) 
Using Eq. ( 12. 91) . we calculate 

(l-p) 4 (2 + 3p-4p 2 -p 4 +p 5 ) 

~ 2(1 - 2p* + 6p3 - 6 p* + 2pt) ■ (5 ' 14) 

For small p, (n)G[K 2 ,jn] has the Taylor series expansion 

(n) G [K 2 ,M = l- h -p + 2p 3 + 7 -p 4 -p 5 + 0(p 6 ) . (5.15) 
For p near to 1, Eq. (I5.14p has the Taylor series expansion, in terms of the expansion variable 

s = l-p, (5.16) 

(n) G[K2M = l -s A + 2s 5 - 2s 7 + 4s 8 + 0(s 9 ) . (5.17) 

For the calculation of the divergence in the correlation length we find, for the ferromagnetic 
Potts model on this strip, 

a 4 

£gik 2 ,M ~ — as a -> oo . (5.18) 
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Hence, setting q = 1 and v = v p , the divergence in the corresponding £ in the bond percola- 
tion problem is 

^G[K 2 ,jn] ~ {1-pY ^ P ^ 1 ' ( 5,19 ' ) 

so that 

vg\k % M = 4 • (5.20) 

For the various infinite-length, finite-width K-regular lattice strips for which we have 
carried out calculations, except for the L y = 2 square-lattice strip with toroidal boundary 
conditions, which involves double bonds in the tranverse direction, we find that the corre- 
lation length in the ferromagnetic Potts model diverges like £ oc a K_1 jq as a — > oo. Thus, 
setting q = 1 and v = v p , this yields, for the corresponding bond percolation problem on 
these strips, £ oc 1/(1 — p) K ~ l as p /* 1. Further calculations with wider strips and other 
lattice types are necessary to determine how general this formula is. The fact that £ diverges 
in the same way for the L y = 2 square-lattice ladder strip and the L y = 2 toroidal square- 
lattice strip indicates that for these strips, doubling the bonds along a direction orthogonal 
to the longitudinal direction does not change the critical behavior of the bond percolation. 
This is understandable, since it is only the longitudinal direction in which the L x — > oo and 
the infnite percolation cluster forms for p = 1. Our general result in Eq. ( I3.25P shows that 
if not only the transverse bonds, but also the longitudinal bonds are doubled, this changes 
v = 2 to v = 4. 

As in our earlier work 17], although in an analysis of physical percolation one is primarily 
interested in real p G [0, 1], it is also useful to investigate the analytic properties of (n) more 
generally in the complex p plane. The reason for this is that singularities in complex p can 
have an important influence on series expansions in small or large p . Here we observe 

that (n) has singularities, which are simple poles, at the zeros of the denominator of Eq. 
f)5.14p . There are five such poles, which we list below to the indicated accuracy: 

{-0.418530 , 0.300885 ± 0.674465? , 1.408380 ± 0.454693i} (5.21) 

Of these, the first is the closest to p = and determines the radius of convergence of the 
small-p Taylor series for (n). 

Thus, studies of percolation quantities on quasi-one- dimensional lattice strips in jl3] and 
here yield valuable insights into the influence of unphysical singularities in series expansions 



about p = and p = 1 for percolation on higher-dimensional lattices [36j. As in [17], we 
can understand these poles in (n) more deeply by noting that although the free energy 
f({G}, q, v) of a given infinite-length (m — > oo) limit of a lattice strip graph only depends on 
a dominant eigenvalue of the relevant transfer matrix, the partition function for a finite strip 
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FIG. 1: Plot of the boundary B qn in the complex p plane for the infinite-length limit G[K2,jn]. 

graph G m is a sum of m'th powers of these eigenvalues. With the substitution v = v p in Eq. 
(12. 7p . these eigenvalues become functions of p. In the infinite-length limit, partition function 
zeros in the complex p plane merge to form boundaries separating regions where a given 
eigenvalue is dominant, i.e. has the largest magnitude and hence determines the resultant 
f({G}, q, v p ). Plots of the resultant boundaries were shown for various infinite-length limits 
of lattice graphs in [13]. For the G[(K 2 ) m , jn, FBC] strip, there are two such A's, given 
above in Eq. (15. lip . Evaluating these for v = v p and q = 1, we have 

^G[K 2 ,jn],+ = JY~p~jE ( 5 - 22 ) 

and 2 

^G[K 2 ,jn],- = p _ p y ■ (5.23) 

The boundary curve is the set of solutions of the equation of degeneracy in magnitude of 
these As, namely 

2|p 2 (l - p) 3 | = 1 . (5.24) 

The solution is a closed egg-shaped curve, shown in Fig. [TJ This curve crosses the real p axis 
at p — —0.418530... (the only real root of the equation 2r 2 (l— r) 3 — 1 = 0) andp = 1.584080... 
(the only real root of the equation 2r 2 (l — r) 3 + 1 = 0). This curve thus constitutes the phase 
boundary in the complex p plane, separating this plane into two regions. As we proved in [13] 
for the infinite-length limit of an arbitrary strip graph, the physical real interval < p < 1 
lies entirely inside the inner region bounded by this curve. The three poles of (n) in Eq. 
( I5.14p lie on this boundary curve. 
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D. (n) for G[K 3 ,jn] 



One can also compute (n) for the m — > oo limit of the family G\{K^) m , jn, BC], viz., 
G[K 3 ,jn}. From our calculation of the partition function for this graph in [34], we obtain 
f(G[K 3 ,jn],q,v) = (1/3) ln[X G [ K3dn ] )0imax \, where X G [K 3 ,jn],o,max is the root of maximal mag- 
nitude, in the relevant interval p G [0,1], of the cubic equation displayed (for the Tutte 



polynomial equivalent to the Potts partition function) in Eqs. (A.6)-(A.9) of 34|. The re- 
sulting expression for (n) is too lengthy to include here, but we will give the resultant Taylor 
series expansion for p near 0, 
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\n) G[ K 3 ,jn] = 1 - 4p + ^-p 3 + 24p 4 + 39p 5 + 0{p 6 ) . (5.25) 



The small-p Taylor series expansions of (^){g> have the general form 13] 



(n) {G} = l-(^)p + .... (5.26) 

where K e ff is the effective vertex degree for {G} and ... denote terms that are higher-order 
in p. Our Taylor series expansions of {n) G \K T ,in\ f° r the r = 2 and r = 3 cases, as well as 
the elementary r = 1 case are in accord with this general form, since K e ff = 3r — 1 for 
these strips, as given in (I5.3p . The increase in K e ff with r means that for a given p, there is 
an increased probability of forming larger clusters, which, in turn, decreases the number of 
clusters per vertex. 

These results on the infinite-length G[K r ,jn] families containing K r subgraphs thus pro- 
vide further insight into percolation on various families of graphs. Parenthetically, we note 
that rather than repeated K r subunits, one could consider bond percolation on a single K N 
graph [l5j]. The graph is much more highly connected than a graph in either of the 
families G[(K r ) m , jn, FBC] or G[(K r ) m , jn, PBC], since for p — 1, each vertex of Km starts 



out connected to every other vertex. Bond percolation on was studied in [37| (reviewed 
in 0]), and it was shown that (in a probabilistic sense) the size of the largest connected 
component on K N diverges as — )■ oo if p > 1/N. 

VI. ANALYSIS OF FAMILIES OF GRAPHS WITH BOUNDED DIAMETER 

A. Motivation 

The essence of bond percolation on a usual lattice A of dimension d > 2 is that as the 
bond occupation probability increases through the critical threshold value, p Ci A, an infinite 
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percolation cluster appears, linking vertices that are arbitrarily far apart on A. Although 
p Ci A = 1 for a one-dimensional or quasi-one- dimensional lattice, the same statement applies. 
One may ask how percolation quantities would behave if one considered families of graphs 
that have the property of a bounded diameter 16] as A" — y oo. For these families of graphs, 
the usual notion of a critical p c beyond which there is a percolation cluster connecting two 
vertices arbitrarily far apart is clearly not applicable. But how would the usual quanties such 
as (n), P(p), and S(p) behave on such families of graphs? Here we address this question 
and calculate exact analytic expressions for these quantities on two families of graphs with 
bounded diameter as A" — y oo. We note that, in addition to the property of bounded 
diameter, both of these families consist of planar graphs which also share a related property, 
namely that they both contain a vertex whose degree goes to infinity as A" — y oo. The second 
family is also self-dual. 



B. Star Graphs 

A star graph Sn consists of one central vertex with degree N — 1 connected by bonds with 
N — 1 outer vertices, each of which has degree 1. (The context will make clear the difference 
between this symbol and the symbol Sn for the symmetric group on N objects.) The graph 
S*2 is degenerate, in the sense that it has no central vertex but instead coincides with L 2 . 
The graph S% is nondegenerate, and coincides with L 3 , while the S n for A^ > 4 are distinct 
graphs not coinciding with those of other families. From the calculation of Z(S n ,q,v,w) in 



24|, we have 

f({S},q,v,w) = M*s) (6.1) 

where 

Xs = q + w — 1 + wv . (6.2) 

Setting v = v p and carrying out the differentiations to calculate (n), P(p), and S(p), we 
obtain 

(n) {s} = l- p , (6.3) 

P(p) {S }=P, (6.4) 

and 

S(p) {s} = l-p. (6.5) 

We will comment on these results after calculating the corresponding quantities for two other 
families of graphs. 
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C. Families of Self-Dual Graphs 



We next consider two families of planar self-dual (SD) graphs. One family is constructed 
by taking a path graph with N — 1 vertices, adding one external vertex, and connecting all 
of the vertices of the path graph to this external vertex with single bonds, except for the 
vertex at one end, which is connected to the external vertex with a double bond. A second 
self-dual family is constructed by taking a circuit graph Cn-i with N — 1 vertices, adding 
one external vertex and connecting all of the N — 1 vertices to this external vertex, thereby 
forming the wheel graph Wh^. In 38j we called these dual boundary conditions (DBC) 
DBC1 and DBC2, and we calculated the Potts model partition functions for them. In the 
iV — > oo limit, these yield the same reduced free energy. From this we computed 17J] 

(1 - p) 3 

( n )(iD) SD = :— 2 • ( 6 - 6 ) 

1 — p + p 2 



(See Fig. 4 of [Tq for a plot.) 

Here we generalize this analysis to the case of a finite external magnetic field in order to 
calculate P(p) and S(p). As was true of the zero-field case, in the relevant limit, namely 
N — > oo, the free energy per vertex is the same for DBC1 and DBC2 self-dual (SD) boundary 
conditions. We give the partition function for the case of self-dual boundary conditions of 
type 2 (DBC2), Z(Wh n+ i, q, v, w), in the appendix. Taking the iV — > oo limit of this family, 
we calculate 

f((lD)sD, q, v, w) = ln[A(i D)s J , (6.7) 

where \(id) sd = ^z,g d ,i,o,i is given as the solution with the + sign in front of the square 
root in Eq. (18. 8p of the appendix. Evaluating Eqs. (" 12 . 1 2[) and (12 . 13[) for this case, we find 

p/ \ p{l-p 2 +p 3 ) 

P (P)(ID) SD = ~tz ~ 2 x 2 (6-8) 

(i — p + p z y 



S(P)(iD) SD = , _ n - o, 3 • (6.9) 



and 

'1 — p) 3 (l +p — p 2 ^ 

[i - p + p 2 y 

These are plotted in Figs. [2] and [3j 

We now comment on the behavior of these quantities for both the iV — > 00 limit of the star 
graph family, {S} and of the self-dual graphs, (1D)sd- I n accord with the general discussion 
given above, (n) is a monotonically decreasing function of p, decreasing from (n) = 1 at 
p = to (n) = at p = 1. Recall that P(p) for the infinite ID lattice vanishes identically 
for p < p c = 1 and has a jump discontinuity to the value P(l) = 1. In contrast, P(p){s} 
and P{p)(iD) S D are both nonzero and monotonically increasing in the interval p G (0,1]. 
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FIG. 2: P(p) for the the N -> oo limit of the self-dual ID graph, {ID) SD . 
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FIG. 3: S{p) for the the N -> oo limit of the self-dual ID graph, (1D) SD . 

Moreover, for the ID lattice, S'(p) diverges as p approaches p c — 1 from below. In contrast, 
decreases monotonically from 1 to as p increases from to 1, while for the same 
range of p, S(p)(iD) SD , starts at 1, first increases, reaches a maximum (of approximately 
1.0545) at an intermediate point (namely, p = 0.103657, a root of the equation dS{p)/dp = 0) 
and then decreases to zero as p /* 1. The differences in behavior with respect to both the 
infinite ID lattice and higher-dimensional lattices show the effect of the fact that the {S} 
and (1D)sd families have bounded diameter. Among these families, one may also discern 
rather different analytic properties in the complex p plane. The various quantities (n), P(p) 
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and S(p) are entire functions for the {S} family. In contrast, (n), P(p), and S(p) for the 
the infinite-length limit (ID) $d have pole singularities at the same two complex-conjugate 
values of p, namely 

p = e ±in/3 , (6.10) 

where 1 — p + p 2 vanishes. The average per-site cluster number (n) has single poles at each 
of these points, while P(p) has double poles, and S(p) has triple poles. In all cases, these 
singularities set the radius of convergence of small-p Taylor series about the origin as unity. 
It is interesting that this is the same as the radius of convergence of the small-p series for 
S(p)iD, but instead of a pole at p — 1, one has here a complex-conjugate pair of poles on 
the unit circle. 



VII. CONCLUSIONS 



In conclusion, in this paper we have presented some exact results on bond percolation. In 
one part of our work we have derived a relation between a percolation quantity on an £-fold 
bond-inflated lattice and the corresponding quantity on the original lattice A evaluated 
with a transformed bond occupation probability pe given by Eq. (13.61) . This is applicable 
for arbitrary lattices and, more generally, A^ — » oo limits of families of graphs, {G}, and 
provides a precise measure of how percolation quantities change as a consequence of this 
bond inflation. As an application of this general relation, we have calculated threshold bond 
percolation probabilities on various bond-inflated lattices. We have shown that this bond 
inflation leaves the universality class of the percolation transition invariant for dimension 
d > 2 but changes it on a one-dimensional lattice and on quasi-one-dimensional infinite- 
length, finite-width strips. This was demonstrated via changes in both the critical exponents 
7 and v. We have also presented expressions for the average cluster number (n) per vertex 
for the bond percolation problem on the infinite- length limits, m — > oo, of a family of highly 
locally connected graphs, namely G[(K r ) m , jn], for several r values. These add to one's 
knowledge of the dependence of percolation quantities on vertex degree and also give further 
insight into singularities of these quantities in the complex p plane. Finally, we have studied 
some families of graphs with bounded diameter and have investigated, via analytic results, 
how this property affects quantities such as (n), P(p), and S(p). 
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VIII. APPENDIX 

Here we give the partition function for the g-state Potts model on a ID graph with self- 
dual boundary conditions of type 2 (DBC2), which is a wheel graph Wh n consisting of n — 1 
vertices forming a circuit graph C n , each connected by a bond to one central vertex. The 
general form of this partition function for a DBC2 strip of length L x = m vertices and width 
L y vertices, denoted Gn, L x x L y (and having n = L x L y + 1 vertices in total) was derived in 



(Eq. (7.17) of) Ref. [22]. It is 



Ly+l n zh (L y ,G D ,d) L y n Z h(L y ,d) 

Z(G D ,L y xL X} q } v,w)=J2^ d) E ^z,G D ,L u Aj) m +™J25 {d) E ( ^ dM T , 

d=l j=l d=0 j=l 

(80) 

where the numbers n Z h(L y ,GD,d) and n Z h(L y ,d) were given for general L y and d in [22| 
and the coefficients k^ and c^ d ' are defined as follows: 

« M = E(-i) i ( 2d " 1 " J ')(g-i) d - i (8.2) 

j=o \ 3 / 

and 

& = £(-lY( 2d - j )(q-l)^ . (8.3) 
j= o V 3 J 

In the general notation of [22], the wheel graph is Wh^ = Gd,L x x L y with L x = N — 1 
and L y = 1. For the family of wheel graphs we thus need the coefficients — q — 1, = 
(q -l){q- 3), c(°) = 1, and = q - 2. We have (from Table 5 of [22]) n Zh (l, G D , 1) = 3, 
n Zh {l,G D ,2) = 1 and (from Table 1 of Q) n zh (l,0) = 2, n Zh (l, 1) = 1. The three 
\z,G D ,i,ij> 3 = 1' 2; 3, are the roots of the following cubic equation: 

£ 3 + a 2 £ 2 + ai C + a = , (8.4) 

where 

a 2 = — v 2 — 3v — q + 1 — u?t> — w (8.5) 
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a i — v(qw + wv 2 + Awv + q + w — 1 + v 2 + qv) 



and 



a = —wv 2 {v + l)(v + q) . 
The two Xz,G D ,i,o,j are the roots of a quadratic equation, 



(8.6) 
(8.7) 



A 



Z,G D ,l,0,j 



2 L 



q+v-l + w(v + iy± {v + q-l + w(v + l) 2 } 2 -Awv (v + l)(v + q) 



1/2 



where j — 1, 2 corresponds to the ± sign. Further, as a special case of the general structural 
results of [22] , 

Xz,g d , 1,2,1 = u , Az,G D ,l,l,l = u • (8.9) 

Thus, explicitly, 

3 2 

q, v, w) = (q-1) J^i^GMjf-'Hq-iyq-^v^+w J^C^G^f^Hq^wv 



N-l 



3=1 



3=1 



(8.10) 

It is easily checked that in the special case of zero field, w = 1, this reproduces our calculation 
of ZiWfiN, q,v) in Js^]. For the present application, we only need the dominant A in the 
region p G [0, 1], i.e., via Eq. ( 12. 71) . v > 0, with ft, — )■ + , i.e., w — )■ 1 + , and this is Az,g d ,i,o,i- 
With the substitution v = v p , w — > 1 + , we denote this simply as X(id) sd in Eq- (16.71) . 

Moreover, on the connection between percolation and the q = 1 Potts model, we note 
that from our exact calculation of the relevant transfer matrix in [22j and the resultant 
Z(L n , q,v,w) and Z(C n ,q,v,w) in [25j], it follows that 



/({£}, 9, w) = f({C}, q, v, w) = ln(A Li i A +) 



(8.11) 



where 



Al,i,o,= 



q — 1 + v + wCl + v) ± 



{q — 1 + v + w(l + v)} 2 — Avw(q + v) 



1/2 



(8.12) 



It is readily verified that applying the differentiations in (12 . 12j) and ( 12 . 1 3[) reproduce the 
known results (13.351) and (13.421) . This derivation is complementary to the usual one via 
combinatoric arguments. 
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